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Abstract 

We demonstrate the existence of a novel superconducting state in high quality two-component 
MgB2 single crystalline superconductors where a unique combination of both type-1 (Ai/^i < 
l/\/2) and type-2 {\2/i2 > l/\/2) superconductor conditions is realized for the two components of 
the order parameter. This condition leads to a vortex- vortex interaction attractive at long distances 
and repulsive at short distances, which stabilizes unconventional stripe- and gossamer-like vortex 
patterns that we have visualized in this type- 1.5 superconductor using Bitter decoration and also 
reproduced in numerical simulations. 
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Vortices in superconductors are characterized by a normal core with size ^ (coherence 
length) and supercurrents flowing over a distance A (penetration depth). If two vortices 
are generated in a type-1 superconductor the normal cores would overlap flrst, due to the 
larger value of C, with respect to A, thus leading to a gain in the condensation energy 
and, consequently, to vortex- vortex attraction l|,l2|. Two vortices in a type-2 material, 
however, would have flrst their supercurrents overlapping, in view of the bigger A, leading 
to vortex-vortex repulsion. An attractive vortex-vortex interaction results in the formation 

n 

of macroscopic normal domains in the intermediate state [3|], while vortex- vortex repulsion 
leads to the appearance of the Abrikosov lattice. {4] 



The recent discovery of the two-component superconductor MgB2|5|], with two weakly 
coupled coexisting order parameters, \&i = "^j^ and \l/2 = "^a, has opened remarkable new 
possibilities both for fundamental research and applications. Among the new research top- 
ics we flnd semi-Meissner stateldl, the violation of the London law and Onsager-Feynman 

Q n 

quantization [7|], non-composite vortices [8|, intrinsic Josephson effect with low-energy in- 
terband Josephson plasmons J9| , mapping physics of two-component superconductors on 
two-condensate Bose systems[10|], superflui dity in liquid metallic hydrogen etc. The 



two-component character of MgB2 
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l4 . [isl is related with two different types of 



electronic bondings, vr and a, giving rise to two different superconducting gaps with en- 



ergies Aj,{0) = 2.2 meV 
the BCS expression ^(0) 
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17j and A^(0) = 7.1 me V 



20[ |. respectively. Using 
hvp/7rA{0), where vf is the Fermi velocity (5.35x10^ m/s 



for the vr-band and 4.40x10^ m/s for the u-band) 21[ for the vr- and cr- bands separately, 
we obtain the two coherence lengths: C,tt{0) = 51 nm and Co-(O) = 13 nm. The calcu- 
lated ^7r(0) value is in agreement with C,n{^) = 49.6 ± 0.9 nm obtained from the flt of 
the vortex proflle measured by scanning tunneling spectroscopy (STS) [l^. The London 
penetration depths A^r and A^r can be found from the respective plasma frequencies UpT, 



and ujpcr 2^: A7r(0) = 33.6 nm and Acr(O) = 47.8 nm. As a result, at least in the clean 
limit, the vr- and a- components of the order parameter in MgB2 are in different regimes: 
= A^(0)/^^(0) = 0.66 ± 0.02 < l/v^=0.707 (type-1) (assuming A^(0) = 33.6 ± 0.9 nm) 
and Ko- = Ao-(0)/,^o-(0) = 3.68 > 1/V2 (type-2). However, when MgB2 is doped with im- 
purities, both components of the order parameter (due to an increase of A and a decrease 
of ^) will fall inevitably into the type-2 regime. Therefore, clean MgB2 represents an excel- 
lent candidate to search for a new type of superconductivity, neither of type-1 nor type-2 
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character, which we coined as type-1.5 superconductivity. In type-1.5 superconductors the 
vortex-vortex interaction is the result of the competition between short-range repulsion and 
long-range attraction and it is expected that vortices could form unusual patterns at low 
applied fields. [6|] 

In this Letter, we present experimental observation of vortex patterns at low vortex 
densities in high quality MgB2 single crystals. The vortex patterns are compared with the 
results of molecular dynamics simulations based on a two-gap Ginzburg-Landau (GL) theory 
which results in peculiar equilibrium vortex structures such as gossamer-like vortex patterns 
and vortex stripes. 

The MgB2 single crystals were grown using a high pressure method as described elsewhere 



231]. The critical temperature of the samples is typically 38.6 K as determined from the ac- 
susceptibility response at zero field. The temperature dependence of the lower and upper 
critical fields Hci(T) and Hc2{T) were determined by means of magnetization measurements 

24, 2^ and by ac-susceptibility measurements, respectively. From the extrapolated value 
Hc2{0) = 5.1 T we obtain ^ab{0) = 8.0 nm close to the value found from BCS theory for .Co-(O). 
In order to estimate the penetration depth of the cr-band we use the theoretical expression 
for the field at which the first vortex will penetrate a two component superconductor [Eq. 



(4) in Ref. 6|]. Considering A7r(0) = 33.6 nm 22| for the type-1 component and taking into 



account C,ab{^) = 'Co-(O) = 8.0 nm and -ffci(O) = 0.241 T from our measurements we obtain 
Aafe = Ao- = 38.2 nm for our samples. 

Coexisting two interpenetrating weakly coupled order parameters is a new quantum state 
with superflow of the two components through each other without any friction. A vortex- 
vortex interaction is derived from the GL theory by numerically minimizing the free energy 



of two vortices with a variational procedure 



261 ]. We use the two-band GL functional: 



FcLl'^a, ^n, A]=F, + F^+ [ d'r[^{^:^^ + *;^,) + -^(V x A) 

J ^ ^fJ'O 



where 



rfV[2E,„]vl/j2 + \E,^\\^^\^ + ^W^l(-zV + ^A)^^\'] (2) 

where a = a, vr and $o is the flux quantum. The estimations for the intrinsic condensa- 
tion energies E^a and the coupling energy E^ are taken from Ref. {2^] and we use the values 
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of and obtained for our samples. The result of the minimization shows that the inter- 
action between vortices is short-range repulsive and weakly long-range attractive, similarly 
to Ref. 6|]. We model a system of overdamped vortices by molecular dynamics simulations. 
The equation of motion for a vortex z is Fj = F™ + ¥j = rjVi, where F^^ represents the 
vortex-vortex interaction and FJ the thermal stochastic force satisfying {Fj(t)) = and 
(Ff(t)Fj(t')) = 27]6ij6{t -t')kBT. r]{=^oHc2/ pn) is the viscosity, Hc2=5.1 T is the upper 
critical field and p„=0.7 fiQ-cm 27| is the normal state resistivity. The system size used 
in our simulations is 2000x2000A^(0). Two systems with number of vortices = 150 and 
400 are initially prepared in a high temperature molten state and then slowly annealed to 
T = 4.2 K with 1000 temperature steps. We let the system stabilize during 2000 time steps 
in each step of temperature. In order to compare the vortex structure of a type- 1.5 with a 
conventional type-2 superconductor we also carry out this simulation for NbSe2 using the 
superconducting parameters from Ref. 



Bitter decoration experiments on MgB2 single crystals were performed at 4.2 K after 
cooling down the sample from above Tc in the presence of an applied field perpendicular to 
the sample surface (field cooling, PC). In this way, a homogeneous vortex distribution all 
over the sample is expected. Details about the technique can be found elsewhere |29|. A 
Bitter decoration image at H = 1 Oe shows clear evidence of an inhomogeneous distribution 
of vortices [Fig.[T](a)] reminding gossamer patterns: local groups of vortices with intervortex 
distances shorter than the average vortex distance (2$o/"\/3-B)^^^ ~ 5/im are separated 
by randomly distributed vortex voids with size of few micrometers. This is in striking 
contrast with the conventional homogenous vortex pattern formed in NbSe2 single crystals 
[see Fig.[T]^b)]. The observed clustering of vortices in MgB2 samples is consistent with 
the theoretical modeling made in Ref. 6| for a two-component superconductor in the semi- 
Meissner state. In Fig.[2]^a) the vortex positions in a selected region of the Bitter decoration 
image shown in Fig.lU^a) are indicated as white dots, while in Fig.[2](b) the results from the 
numerical simulations for a two-component superconductor are shown. We also calculate 
the vortex configuration for a reference conventional type-2 superconductor. The obtained 
vortex structure, considering A = 69 nm and ^ = 7.7 nm 28|, is similar to the one observed 
in NbSe2 samples [see Figs.[2](c) and (d)]. 

In order to characterize the inhomogeneity of the vortex structure (VS) we calculate the 
distribution of first neighbor distance. Pa [Figs.[2]^e) and (f)]. The first neighbor distance. 
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a, is calculated by means of the Delaunay triangulation of the vortex structure. For the 
VS of NbSe2, the distribution is Gaussian with a relative standard deviation 6 = SD/ < 
a >= 0.224, where SD is the standard deviation of the Gaussian fit to the experimental 
data and < a > is the average first neighbor distance. On the other hand, in MgB2 samples 
Pa is quite broad as a consequence of the inhomogeneous arrangement of vortices and has 
additional peaks at distances shorter and longer than the most probable separation [see the 
red and green arrows in Fig.[2]J^e)]. The distribution of the first neighbor separation of the 
vortex structure obtained by simulations of a two-component material has also a three-peak 
structure which is fairly broader when compared to the one obtained in the case of a one- 
component conventional type-2 superconductor [Fig.[2]^f)]. The peaks at short distances in 
Figs.[2](e) and (f) correspond to an average minimum separation between vortices. 

Figure[3t^a) shows a Bitter decoration image of the VS in MgB2 crystals at = 5 Oe and 
T = 4.2 K. The vortex distribution appears also to be inhomogeneous at this field but in 
a rather different manner than in the experiment described above. In some regions of the 
sample, vortices agglomerate forming stripes while voids are formed on considerably larger 



cirGcLS. 



30l | The present experimental results seem to be in contradiction with previous reports 



of the vortex structure at low fields in MgB2 samples |3l)]. However, it is important to note 
that in Ref.lSlI the authors show an image of the VS at 4 0e in a very small region of the 
sample (approximately 10 x 10 yum^). Therefore, it is not possible to determine whether the 
vortex distribution in the samples studied in Ref. 3l| is uniform all over the sample or not. 

Although the vortex stripe pattern is rather disordered it is still possible to determine 
an average direction of vortex stripes as the one defined by the dashed yellow lines in 
Fig.[3](a). Considering this definition, we calculate the vortex density in lines parallel to 
the vortex stripes as a function of the distance measured along the direction of the yellow 
arrow in Fig.[3]^a). In the inset of Fig.[3](b) we plot the linear vortex density normalized 
by the average value for both the MgB2 and NbSe2 VS at 5 Oe. In the case of MgB2, 
fluctuations of the vortex density of the order of 50% are observed. A similar calculation 
along lines perpendicular to the stripes shows that the standard deviation of the mean 
value is of the order of 30%. The large fluctuations of the vortex density in MgB2 are in 
contrast to what is observed in NbSe2 crystals where the standard deviation of the vortex 
density is approximately 1% of the average value. A remarkable similarity is found between 
the experiments and simulations at still low density but higher than the one shown in 
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Fig. [21 Disordered vortex stripes are formed in the two-component superconductor while a 
homogeneous distribution is apparent in the case of a conventional type-2 material [Figs.[3](c) 
and (d), respectively]. Consistently, the vortex density is seen to fluctuate in the direction 
perpendicular to the vortex stripes in the two-component material, as shown in the inset of 

Fig.EKd). 

Composition analysis (via an electron microprobe in a field emission scanning electron 
microscope) in an area across the stripes [in the direction of the white arrow in Fig. 3(a)] 
shows no significant variations in Mg or B content, thus ruling out the possibility to at- 
tribute the stripe formation to inhomogeneous surface pinning distribution. There is also 
no observed correlation between the vortex positions and localization of microdefects. 

Bitter decoration images at H = 0.5 Oe also show a very inhomogeneous vortex distribu- 
tion (not shown). On the other hand, at = 10 Oe the VS in MgB2 samples is similar to the 
one in NbSe2 crystals indicating that this novel superconducting phase in two-component 
y^pe-l.S superconductors is only accessible at very low applied fields as predicted in Ref. 

i 

In conclusion, the type-1.5 superconductivity is a totally new state which combines two 
regimes (type-1 and type-2) in the same sin gle m aterial (clean MgB2 and possibly also other 
two-gap materials such as Bao.6Ko.4Fe2As2 3^)- The vortex matter in type-1.5 supercon- 
ductors behaves in an extremely unusual way. The combination of the vortex-vortex repul- 
sion and attraction in the same material leads to the appearance of novel vortex patterns: 
gossamer-like vortex arrays and vortex stripes. Both novel patterns have been directly visu- 
alized by Bitter decorations on high quality single crystals. Moreover, analytical modeling 
of exotic vortex-vortex interaction (attractive mixed with repulsive) and extensive molecular 
dynamic simulations are in good agreement with our experimental data. 
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FIG. 1: Magnetic decoration images of the vortex structure at T = 4.2 K and H = 1 Oe. The 
vortices in the MgB2 single crystal (a) are distributed inhomogeneously while in the NbSe2 single 
crystal (b) they are arranged forming a disordered Abrikosov lattice. The scale bars in the images 
correspond to 10 /im. Notice that the density of vortices in the decoration experiments represents 
the internal field B rather than the applied field H. This leads to different number of decorated 
vortices for NbSe2 and MgB2 , even at the same applied field. 

FIG. 2: (a) Experimental vortex locations in a selected part of the image shown in Fig. 1(a). The 
vortex configuration resulting from the numerical simulations in a two-component superconductor 
at low density is shown in (b) evidencing an inhomogeneous spatial distribution of vortices. In 
both cases, the regions enclosed by the dashed white line indicate voids of vortices caused by the 
inhomogeneous distribution. In (c) the vortex pattern obtained by a magnetic decoration of the 
NbSe2 crystal at 1 Oe is shown and (d) corresponds to the vortex pattern obtained by a numerical 
simulation of a type-2 superconductor. The white scale bars correspond to 10 ^m. (e) and (f) 
display the distribution of first neighbor distance, Pa, of the experimental and theoretical vortex 
structures, respectively. The distributions obtained for the simulations of a type-2 material and 
for the NbSe2 vortex structure are Gaussian. In the case of MgB2 Pa shows additional peaks at 
distances shorter and longer than the most probable separation (see the red and green arrows). 
Pair of vortices separated at the distances where the additional peaks are located are highlighted in 
(a) and (b) by red and green circles. The light blue circles correspond to pair of vortices separated 
by the most probable separation. 
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FIG. 3: (a) Magnetic decoration image showing a stripe-like vortex pattern in the MgB2 single 
crystal at H = 5 Oe. (b)Disordered Abrikosov lattice with a homogeneous vortex density obtained 
at -ff = 5 Oe in the NbSe2 sample. The formation of vortex stripes is also observed by numerical 
simulations of a two-component type-1.5 superconductor (c) in contrast to a homogeneous vortex 
distribution in a type-2 superconductor at the same vortex density (d). The scale bars in the 
images correspond to 10 fim. Vortex density along lines parallel to the vortex stripe direction 
(yellow dashed lines in (a)) for MgB2 and NbSe2 vortex structures. The variation of the vortex 
density as a function of the distance measured along the direction perpendicular to the stripes 
(yellow arrows in (a) and (b)) is shown in the insert of (b). The curves are normalized by their 
respective average density. The results of a similar calculation performed on the simulated vortex 
structures are shown in the inset of (d). 
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